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Comment on " Regional Versus Global Entan- 
glement in Resonating- Valence-Bond States" 

In a recent Letter Chandran and coworkers study 
the entanglement properties of valence bond (VB) states. 
Their main result is that VB states do not contain (or 
only an insignificant amount of) two-site entanglement, 
whereas they possess multi-body entanglement. Two ex- 
amples ("RVB gas and liquid") are given to illustrate 
this claim, which essentially comes from a lower bound 
derived for spin correlators in VB states. While we do 
not question that two-site entanglement is generically 
"small" for isotropic VB states, we show in this Com- 
ment that (i) for the "RVB liquid" on the square lattice, 
the calculations and conclusions of Ref. [l[ are incorrect, 
(ii) A simple analytical calculation gives the exact value 
of the correlator for the "RVB gas", showing that the 
bound found in Ref[l| is tight. ( iii) The lower bound for 
spin correlators in VB states is equivalent to a celebrated 
result of Anderson dating from more than 50 years ago. 

The SU (2) symmetry of VB states guarantees that any 
two-spin reduced density matrix is a "Werner state" fully 
characterized by a parameter p. The considered pair of 
spins is entangled if p > 1/3. Chandran et al. used 
quantum information concepts such as monogamy of en- 
tanglement and quantum telecloning to obtain bounds 
on p. The number p is simply related to the correla- 
tor (Si.Sj) between these two spins 1/2 (S = tr/2, with 
a Pauli matrices). We have {Si.Sj) = —3/4p (and not 
"exactly equal to the parameter p" as stated in Ref. [l[). 

(i) The "RVB liquid" is the equal amplitude superpo- 
sition of all nearest-neighbour (NN) VB coverings of a 
bipartite lattice. Exact results can be obtained for small 
sizes L of the square L x L lattice. For L = 4, we do 
not recover the value p ~ 0.2004 of Ref. [l[, but find 
p = 0.4457579115872 for periodic boundary conditions 
(BC) and p — 0.2281115037 in the interior of a sample 
with open BC. However, what really matters is the be- 
havior for large L. Exact calculations are difficult in this 
case, but Monte Carlo calculations are possible We 
computed the NN correlator (Si.Sj) for large samples (up 
to L = 128) on the square lattice with periodic BC. The 
data of Fig. 1 shows that p is larger than 1/3 in the ther- 
modynamic limit (we find p = 0.3946(3), resulting in an 
entanglement of formation of ~ 0.0215). Therefore, the 
"RVB liquid" on the square lattice does possess two-site 
(NN) entanglement, contrary to the claim of Ref. [l|. 

(ii) The "RVB gas" is the equal amplitude superpo- 
sition of all bipartite VB coverings of a bipartite lattice. 
This is in fact the projection into the singlet sector of the 
(magnetically ordered) Neel state on this lattice. This 
observation can be used to calculate p exactly. The total 
spins Sa and Sb on sublattices A and B are maximal, 
couple antiferromagnetically and form a singlet (total 
spin S = 0). For a system of 2N spins, S A = S B = N/2. 


One then easily obtains that (Si.Sj) = -1/4 - 1/(2 JV) 
if i and j belong to different sublattices. The equivalent 
exact result p = 1/3 + 2/(3N) shows that the telecloning 
bound p < 1/3 + 2/(3A^) is tight. Two-site entanglement 
is therefore present in any finite "RVB gas" and vanishes 
only in the thermodynamic limit. 


FIG. 1: Werner parameter p as a function of inverse linear 
system size 1/L for the square lattice "RVB liquid" . 

(iii) The telecloning bound on v in Ref. [l| reproduces 
an inequality found by Anderson [3|] , who derived a lower 
bound for the energy of antiferromagnetic spin models. 
Take a spin at site i, separated by any distance from 
a number z of symmetry-equivalent spins j: (Si.Sj) (as 
well as p) is identical for all z spins at sites j. In this case, 
the telecloning bound is p < 1/3 + 2/(3z) or equivalently 
(Si.Sj) > — 1/4— l/(2z), the result derived by Anderson. 
His result (of variational nature) on correlators in a given 
state is very general: it holds also for states other than 
singlets, is independent of any Hamiltonian and can be 
refined further (see e.g. Ref. [J]). 

In conclusion, the bound obtained with quantum infor- 
mation techniques [l[ has been familiar in the condensed 
matter context for a long time. Nevertheless, it is inter- 
esting to see that it can be derived in a totally different 
framework. For the two examples chosen in Ref. [H, typ- 
ical condensed matter methods allowed us to provide in 
one case an exact solution, and to show that the results 
of Ref. [l| are incorrect in the other one. 

We thank M. Mambrini for useful discussions, and D. 
Kaszlikowski and V. Vedral for correspondence. 
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